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Abstract. The periodic Lorentz gas with external field and iso-kinetic thermostat is 
equivalent, by conformal transformation, to a billiard with expanding phase-space and 
slightly distorted scatterers, for which the trajectories are straight lines. A further time 
rescaling allows to keep the speed constant in that new geometry. In the hyperbolic 
regime, the stationary state of this billiard is characterized by a phase-space contraction 
rate, equal to that of the iso-kinetic Lorentz gas. In contrast to the iso-kinetic Lorentz 
gas where phase-space contraction occurs in the bulk, the phase-space contraction rate 
here takes place at the periodic boundaries. 



Submitted to: J. Stat. Mech. 

PACS numbers: 05.45.-a,05.70.Ln,05.60.-k 

E-mail: fbarra@dfi.uchile.cl, thomas.gilbert@ulb.ac.be 



Letter to the Editor 



2 



Over the last twenty years, the study of time-reversible dissipative systems with 
chaotic dynamics as mechanical models of non-equilibrium stationary processes has 
played an important role in shaping our understanding of the connections between 
irreversible macroscopic processes and the reversible dynamics that underly them at 
the microscopic level, see [U and refs. therein]. Among other large deviation relations, 
the Fluctuation Theorem [21 [3] is a central result of this approach. Related results have 
appeared since then [H |5], which have been put to the test experimentally in small 
circuits and biological molecules, see [6j and refs. therein]. 

The forced periodic Lorentz gas with Gaussian iso-kinetic thermostatting was first 
proposed by Moran and Hoover [7] as such a time-reversible dissipative chaotic model 
for conduction, where irreversibility manifests itself in the fractality of phase-space 
distributions. 

A rigorous analysis of that system was later provided by Chernov et al. [8], 
establishing a relation between the sum of the Lyapunov exponents and the entropy 
production rate of this system. The thermostat is here a mechanical constraint, 
chosen according to Gauss' principle of least constraint, which acts so as to remove 
the energy input from an external field [9]. Under this constraint, the kinetic energy 
remains constant and thus fixes the temperature of the system; no interaction with 
a hypothetical environment is needed in order to achieve thermalization. Rather the 
Gaussian thermostat causes dissipation in the bulk. Due to the presence of a regular 
array of circular scatterers on which non-interacting particles collide elastically, this 
system sustains a chaotic regime, at least so long as the external forcing is not too 
strong. As shown in [8], this model has a unique natural invariant measure, with one 
positive and one negative Lyapunov exponents, whose sum is negative -a signature of 
the fractality of the invariant measure- and identified as minus the entropy production 
rate [10] . The comparison with the corresponding phenomenological expression provides 
a relation between the phase-space contraction rate and conductivity. 

As will be reviewed shortly, the trajectories of the iso-kinetic billiard are integrable 
from one collision to the next. It was shown by Wojtkowski [TT] that there exists 
a conformal transformation on the Complex plane that transforms those curved 
trajectories into straight lines. Formally, this amounts to introducing a metric, specified 
according to the conformal transformation, which yields an identification between the 
Gaussian iso-kinetic trajectories and the geodesies of a torsion free connection, called 
the Weyl connection, for which the metric is preserved under parallel transport. This 
is a natural generalization of the symplectic formalism for the Gaussian iso-kinetic 
motion introduced by Dettmann and Morriss [12] . The corresponding billiard, whose 
trajectories follow the geodesies of the Weyl connection, and undergo elastic collisions 
when they reach the boundary, is referred to as the billiard W-flow. For our purpose, it 
will be sufficient to focus our attention on the conformal transformation itself, leaving 
aside the formal aspects of this geometric construction, which we will henceforth loosely 
refer to as Weyl geometry. 

In what follows, we will present the details of this conformal transformation and 
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study the dynamics of the billiard. Our main observation is that the bulk dissipation of 
the iso-kinetic Lorentz billiard, which accounts for the positive entropy production rate, 
disappears under the conformal map, after time reparametrization. The phase-space 
contraction rate and its identification with the entropy production rate are nevertheless 
recovered because of the periodic boundary conditions, which, under the conformal map, 
induce a phase-space contraction rate which, in average, is equal to the bulk dissipation 
of the iso-kinetic Lorentz billiard. We mention that a billiard with a similar mechanism 
of contraction of phase-space volumes at the borders was considered in [13]. However 
the connection to iso-kinetic dynamics was not discussed there. 

We consider the two-dimensional iso-kinetic periodic Lorentz channel with constant 
external field of magnitude e. Due to the periodicity of this system, the dynamics can be 
studied in a unit cell with periodic boundary conditions, as displayed on the left panel 
of Fig. [U This cell is an open rectangular domain, centered at the origin, with one disk 
at the cell's center, and four others at its summits (±1/2, ±v3/2). All the disks have 
identical radii, which we denote by a, chosen so as to satisfy the finite horizon condition, 
i. e. a/3/4 < a < 1/2. The width of the cell is here taken to be unity, and its height \/3. 
Periodic boundary conditions apply at x = ±1/2 and reflection at y = ±v^3/2. The 
direction of the external field is taken towards the positive x direction. Thus particles 
typically tend to move in that direction, winding around the cell from one boundary to 
the other. 




Figure 1. (left) Unit cell with a typical trajectory, corresponding to the external field 
e = 1; (right) The image of that cell and trajectory under the conformal map Q. The 
initial condition is marked by a circled dot. 

Let (x t ,y t ) denote the coordinates of the particle at time t, and let 6 t denote the 
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velocity angle, measured with respect to the x-axis. Assume p 2 = x\ + y 2 = 1, which we 
take as the definition of the temperature. In between collisions the equations of motion 
are specified by 

X t = COS0 t , (1) 

y t = sm6 t , (2) 
9 t = -esin0 t , (3) 
with the solutions 

1, sinfl 

x t = x + -log^— -, (4) 
e sm 9 t 

yt = yo~\ — , (5) 



6 t = 2 arctan 



tan — exp(— et) 



(6) 



When the trajectory collides with a disk, 9 t changes according to the rules of elastic 
collisions. 

Let z = x + ly and f(z) = —e(x + zy). That is, / is a holomorphic function whose 
real part is specified by the potential associated to the external field. According to 
Wojtkowski the conformal mapping 

F{z) = J VL V [-f{z)]dz, (7) 

takes the trajectories (jHEJ) into straight lines. This is to say F{t) = F(z t ) is a straight 
line in the complex plane. Furthermore, since the transformation is conformal, the 
elastic collisions are mapped into elastic collisions. The periodic boundary conditions 
are modified under this map, as will be discussed shortly. 

The map (jTj) takes a point with coordinates (x, y) to the point (u, v) in the complex 
plane, 

x -\- ly i — > u ~\~ iv — F(x + ly), 

= -{exp[e(x + zy)]-l}. (8) 

Note that the term — 1/e was introduced so as to keep the origin fixed under the 
transformation. 

We show that trajectories are straight lines in the (u, v) plane and identify their 
slopes. To that end, we compute the time derivative of F(zt) with z t = Xt + iyt- From 
Eqs. (HI, 

4^F(z t ) = exp[e(x t + iy t )]^r{x t + iy t ), 
= exp[e(x t +iy t )\ exp(i9 t ), 

= exp(ext) exp[i(9 t + ey t )]. (9) 

Notice that <p t = 9 t + ey t is a conserved quantity, as immediately seen from Eq. (jSJ) . If 
we further change the time variable t to s, such that 

ds = exp(ex t )dt, (10) 
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we can write 

-^-(u s + iv s ) = exp(z0 o ), (11) 

which is to say (u s ,v s ) is a straight line trajectory of slope tan(0 o ), as announced, and 
constant speed. This result may indeed have been anticipated, since the transformation 
to the new set of variables coincides with the one for which the iso-kinetic trajectories 
are transformed into geodesies of the metric exp(2ex)(dx 2 + dy 2 ), see [T2] . 

Thus, in terms of the (w,t>,0) variables, the time-evolution of trajectories between 
collisions is specified according to : 

u s + iv s = u + iv + e^°s (12) 

4>s = 0o- (13) 

The time reparametrization Eq. ( flOj) introduces a new time scale, which is natural 
in the new geometry. However it depends on the details of the trajectory. For future 
reference, we will make a distinction between this natural time s and the physical time 
t. We integrate Eq. ( flOl) in order to express the natural time variable, s, in terms of the 
physical one, t : 

s — [( u o + l/ e ) cos 0o + sin 0o] [cosh(et) — 1] 



+ y> + l/e) 2 + ^sinh(et). (14) 

As opposed to the Gaussian iso-kinetic trajectories, the phase-space volumes are 
obviously preserved along the trajectories Eqs. (11211131) . The time variables are however 
different and, as is clear from Eq. ( TH1) . phase-space volumes are not preserved when the 
trajectories are parameterized according to the physical time t. 

To make the comparison clearer, we first consider the Gaussian iso-kinetic dynamics, 
which maps trajectories (xo,yo,9 ) \— > (x t ,yt,O t ) according to Eqs. (JM6]). This map 
contracts phase-space volumes, with a rate given by minus the logarithm of the Jacobian 
of the application, 

= (cosh tt + cos 9o sinh et) _1 , 

= exp[e(x - ar t )], (15) 

i. e. e(x t — xo) is the phase-space contraction rate of the trajectory taking (xo, yo, #o) to 
(x t ,y t ,6 t ). 

In the Weyl geometry of the billiard, on the other hand, the phase-space contraction 
rate of the physical time trajectory (wo,fo,0o) l— * ( u t, v t,(fit) is minus the logarithm of 
the Jacobian, 



- 1/e) cos 0o + ^os: 
v/K) + l/e) 2 + i 



(m + 1/e) cos 0o + ^o sin 0o . , . 
cosh et H , sinh er, 



V 



v/(u + l/e) 2 + ^' 



(16) 
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This is the inverse of Eq. ( I15p . as easily checked. Therefore, so long as the time scales are 
the same, the phase-space volume contraction changes sign, going from the dynamics of 
the iso-kinetic billiard to that of the W-flow. Moreover minus the logarithm of the latter 
expression yields e(«o — u t) only to first order in e. The reason for this is connected to 
the transformation law of x, Eq. (EJ) : x — 1/elog yju~+ 1/e) 2 + v 2 . 

As remarked by Dettmann and Morriss [12] , the periodic boundary conditions break 
the Hamiltonian structure of the billiard in the Weyl geometry. Indeed, they do not 
preserve phase-space volumes as we now show. 

In order to analyze the billiard trajectories in the Weyl geometry, we must first 
consider the transformation of the unit cell under the conformal map, Eq. ([8]), as shown 
on the right panel of Fig. HJ It is straightforward to check the rectangular cell is mapped 
to a trapezoidal figure with curved lateral sides, and the disks to slightly distorted, 
flattened ones. Here are the parametric equations of these elements : 

upper/lower sides : jw + iv = e _1 [exp(et ± iV3e/2) — 1], — 1/2 < t < 1/2 j, 
right/left sides : ju + iv = e -1 [exp(±e/2 + iVZet) - 1], -1/2 < t < 1/2 j, 
central disk : |m + iv = e _1 {exp[ecr exp(z</>)] — 1}, < <p < 27r|, 

outer disks : jw + iv = e _1 {exp [e( ± 1/2 ± iV3/2 + crexp(«0))] — 1}, 

< 6 < 2ir\. 



Notice that the upper and lower borders are oblique lines of slopes ±V3e/2, 
while the left and right borders are concentric arc-circles of center (— l/e,0) and 
respective radii exp(±e/2)/e. Thus, the periodic boundary conditions induce phase- 
space expansion/contraction as the variables u and v are rescaled by a factor exp(±e) 
when the trajectory crosses from left to right or right to left. 

In order to analyze the ergodic properties of the W-flow, we turn to the definition 
of the Birkhoff coordinates of the Weyl billiard in the periodic cell, which specify the 
evolution of trajectories from one collision to the next (including collisions with the 
boundaries). Three coordinates are necessary here : which denotes the arc-length 
along the boundaries of the unit cell (walls and disks), u>, the modulus of the velocity 
vector, and m, by which we denote the sinus of the angle between the outgoing trajectory 
and the normal to the obstacles' boundary, that is the tangent component of the unit 
vector in the direction of the velocity. As a result, in the stationary state, so long as 
the regime is hyperbolic, there are three Lyapunov exponents which characterize the 
chaoticity of the Birkhoff map. We denote them by Ai > A2 > A3, where Ai and A3 
are associated to the (q, ^-dynamics, and A2 to the cj-dynamics, corresponding to the 
direction of the flow, and which is zero as it will turn out. 

We first consider the dynamics of u. Here we have to make a distinction between 
the two time parameterizations. If, on the one hand, we integrate trajectories with 
respect to the physical time t, the modulus of the velocity is typically expanded in the 
bulk according to Eq. (fl6l) . while it is contracted at the opposite rate by the periodic 
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boundary conditions. That the particle's velocity must be rescaled by a factor e ±e at the 
periodic boundaries can be seen from the expression of the velocity, given by Eq. (fTU|) . 
here expressed in terms of the variables u, v, <p, 

uj t = -j-s = e j [(uq + 1/e) cos O + ^ o sin O ] sinh(et) 

+ ^/( U o + l/e) 2 + ^o 2 cosh(et)}. (17) 

This rescaling yields a contribution to A2 equal tdj] — eJ c , where J c is the (signed) 
average number of boundary crossings per collision, i. e. the average drift velocity per 
collision, a positive quantity. The other contribution occurs from the bulk, because of 
the dissipation that takes place along the trajectories, and is given by the logarithm of 
Eq. (Tl6l) . As we noted by comparing it to Eq. (TT5T) . this contribution is opposite to the 
phase-space contraction rate of the iso-kinetic billiard along the image trajectory. The 
latter quantity is equal to the entropy production rate of the iso-kinetic Lorentz gas, 
given by the product of the applied force by the resulting drift current [8], i. e. eJ c . 
Thus the bulk contribution to A 2 is identical but opposite to the contribution due to 
the boundary conditions. Therefore, 

A 2 = 0. (18) 

If, on the other hand, we integrate the trajectories with respect to the natural time, 
i. e. with respect to s, the dynamics of u are trivial. Indeed the velocity is constant in 
the bulk and remains so under the periodic boundary conditions. This is because the 
s increments are rescaled at the boundaries identically to the length increments, as can 
be seen from Eq. (1171) . Thus, trivially, we again retrieve A 2 = 0. 

As noticed earlier, Eqs. (j 12111 3D . the dynamics of q and w preserves the phase-space 
volumes in the bulk. This remains so whether the time variable is s or t. Since the 
trajectories follow straight lines, the evolution of these two variables from one collision 
event to the next is the same, whether the velocity changes in the bulk or remains 
constant. However, as noted earlier, phase-space volumes are not preserved when the 
periodic boundary conditions apply. The arc-length coordinate is contracted by a factor 
e ±e . We can therefore write minus the sum of the two Lyapunov exponents associated 
to the dynamics of q and w as 

- Ai - A 3 = eJ c , (19) 

The total phase-space contraction rate is the sum of Eqs. (}T8"]) and (JT9~j) . equal to 
the product of the current multiplied by the amplitude of the external field. This is 
obviously the same result as obtained for the iso-kinetic Lorentz billiard [8], which we 
did expect since the winding number, and hence the current J c , are invariant under 
the conformal map. Table [1] summarizes our findings. Notice that the current J c is 
here measured with respect to the collision dynamics. The identification between the 
phase-space contraction rate and the entropy production per unit time (in the regime 

% For definiteness we assume the particle density to be normalized to unity. 
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Dissipation 


Bulk 


Boundaries 


Total 


GIK LG 


eJ c 





eJ c 


W-flow (t) 


-eJ c 


2eJ c 


eJ c 


W-flow (s) 





eJ c 


eJ c 



Table 1. Comparison of the different contributions to the phase-space contraction 
rate for the Gaussian iso-kinetic Lorentz gas and its two equivalent representations 
in the Weyl geometry, whether the integration variable is t, the physical time, or s, 
the natural time, Eq. (fT¥]) . In the former case, the contraction of the velocity at the 
periodic boundaries compensates the expansion of phase-space volumes in the bulk. 

of linear transport) can be done by substituting the average current per collision, J c , 
by the corresponding average current per unit time, which we denote by J t . This is 
straightforward since the two quantities differ only by a conversion factor, the average 
free flight time, or average time between collisions. On the other hand, it should be 
mentioned the conversion to a current measured in the units of s is not meaningful since 
this integration variable changes from one trajectory to another. Each trajectory travels 
with its own clock, which impedes a statistical interpretation. 

To conclude, we have shown how the transformation of the Gaussian iso-kinetic 
Lorentz channel to a Weyl geometry maps the trajectories into straight lines with 
constant speed, so long as we rescale the time variable according to Eq. (flOl) . The 
geometry of the cell is modified under this transformation. The vertical walls of the 
Lorentz channel are here replaced by concentric circles whose radii differ by a factor 
given by the exponential of the external field. The disks are also deformed and can even 
have concavities when the field strength becomes too large (e > 1/er, which corresponds 
to the transition to non-hyperbolic regime, as shown in 

The average phase-space contraction rate, equal (at least in the small e regime) to 
the entropy production rate, is here due to the application of the periodic boundary 
conditions. A remarkable consequence of this, is that the calculation of the phase-space 
contraction, as well as its identification to the entropy production rate, are immediate. 

It should be pointed out that the expression of the phase-space contraction rate, 
Eq. (1191) . makes no reference to the actual time scale since it relies solely on the collision 
dynamics. It is in particular transparent to the change of time scales, Eq. ( flOl) . However, 
as we argued, only the physical time is relevant to statistical averages. In this system 
the velocities change with the physical time and get rescaled at the periodic boundaries. 
The velocity dynamics of the natural time are on the other hand trivial. Thus one 
might wonder what happens when s is taken as the physical time, in which case we 
lose the connection with the Gaussian iso-kinetic dynamics. This is what we do in [14J, 
where we discuss an expanding billiard model similar to this one, but for which the 
physical time dynamics is trivial between collisions and preserves phase-space volumes 
in the bulk. In this case, the speed variable must be rescaled at the periodic boundaries, 
and is associated to a third Lyapunov exponent, here negative. The mobility changes 
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by a factor of two, and one nevertheless retrieves an identification between entropy 
production and phase-space contraction rates in the linear transport regime. 
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